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Abstract. We consider one-dimensional tubes containing bosonic polar molecules. 
The long-range dipole-dipole interactions act both within a single tube and between 
different tubes. We consider arbitrary values of the externally aligned dipolc moments 
with respect to the symmetry axis of the tubes. The few-body structures in this 
geometry arc determined as function of polarization angles and dipole strength by 
using both essentially exact stochastic methods and the harmonic approximation. The 
main focus is on the three-, four-, and five-body problems in two or more tubes. Our 
results indicate that in the weakly-coupled limit the inter-tube interaction is similar 
to a zero-range term with a suitable rescaled strength. This allows us to address 
the corresponding many-body physics of the system by constructing a model where 
bound chains with one molecule in each tube are the effective degrees of freedom. This 
model can be mapped onto one-dimensional Hamiltonians for which exact solutions 
are known. 
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1. Introduction 

Heteronuclear molecules with strong electric dipole moments are a current pursuit of 
the cold atom community [H [21 [3J El 13 E|. The dipole-dipole forces in such systems 
can be externally controlled and give access to long-range and anisotropic interactions. 
Strong interactions can lead to rapid chemical reaction loss which can, however, be 
suppressed by considering low- dimensional trapping potentials [7J [8]. In the case of 
one or several one- dimensional tubes holding dipolar particles theoretical works indicate 
that a number of interesting physical states can be realized. Non-trivial Luttinger liquid 
states P, HOI EH H21 H3] , superfluids, supersolids and stripes pH [15j [16] , liquids of trimers 
and crystals of different complexes [T7J [18], Mott insulators [191 120], exotic quantum 
criticality [211 122], zig-zag transitions [23] and few-body states of several molecules 
pHl 1251 [26] have been discussed in the literature. 

In the current paper we will address the few-body structures that can be expected 
in an array of several tubes. A schematic setup for three tubes is shown in figure [U The 
orientation of the dipolar molecules in the tubes can be controlled by an external electric 
or magnetic field depending on whether magnetic atoms or heteronuclear molecules are 
used. This opens up a host of interesting phenomena since the potentials of two dipoles in 
a single tube and the potential between two dipoles in different tubes change magnitude 
and sign as one changes the angles shown in figure [U 

Here we are concerned with the important question of existence of iV-body bound 
states in multi-tube geometries. In order to address this we need to carefully consider the 
limits of small dipole moment and we consequently develop a perturbative approach to 
weakly bound states in one dimension that can handle the dipolar potential. At stronger 
coupling we consider a harmonic approximation. Both are compared to numerical 
results from the stochastic variational method that allow us to test the accuracy of the 
analytical formalism. A comparison is also made to exact results for iV-body bosonic 
systems interacting through zero-range attractive forces, and we ask whether dipolar 
systems in one dimension approach the exact results in given limits and can in turn 
be effectively described by (properly renormalized) zero-range interactions. This is 
particularly important for many-body studies since zero-range interactions are not only 
convenient to work with but also often allow analytical results to be obtained. We also 
consider the two-body scattering dynamics of the dipoles and compare to zero-range 
results. Finally, we discuss the impact of our results for the many-body physics of 
dipolar bosons in one- dimensional geometries. 

2. Basic setup 

We consider the setup depicted schematically in figure CD, i.e. an array of equidistant 
one-dimensional tubes containing dipolar particles with dipole moments aligned by an 
external field that does not interfere with the tubular geometry. The dipolar particles 
are identical bosons with mass m and dipole moment D. The potential between two 




Figure 1. Schematic of the setup for the case of three equidistant one-dimensional 
tubes along the x-direction in the ccy-plane with a distance between adjacent tubes of 
d. The dipolar moment, _D, points along the direction specified by the two angles </> 
and 6 as defined on the figure. Dark filled circles indicate the dipolar atoms in the 
tubes. The tubes have a thickness given by I. In typical experimental setups the tubes 
also have an in-tube confinement. This is indicated on the figure by the shrinking of 
the tubes at both ends. This in-tube confinement will, however, be neglected in this 
work. 

dipoles, V(n,x), can be written in the form 

md\ lt . rr J + n 2 -3cos(fl) 2 [f cos(0) + nsin(0)] 2 

—-V(n,x) = U- 2 — ; ) (1) 

ft 2 K ' (g +rl 2)5/2 v > 

where d is the distance between two adjacent tubes and n is an integer such that nd is 
the intertube distance (n = 1 for nearest neighbour (adjacent) tubes, n = 2 for next- 
nearest neighbours and so on). Here we have defined the dimensionless dipole coupling 
strength, U = j*nr, where D is the dipole moment (absolute value of the vector D in 
figured]). For n = we get the intratube interaction in the limit where the tubes are 
strictly one-dimensional. In experimental setups, arrays of one-dimensional tubes are 
constructed by applying optical lattices to the dipolar gas [3 [8]. In this case the tubes 
are not strictly one- dimensional but will have some width along the transverse direction 
that is determined by the laser intensity. This can be translated into a gaussian wave 
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packet in the transverse direction that will be increasingly localized in space as the laser 
intensity increases. An effective interaction that takes this into account by integrating 
out a gaussian wave packet can then be obtained [25], which will in the limit of zero 
gaussian width reduce to the expression in (JTJ). We will assume that the lattice is very 
strong so that the strict one-dimensional expression above is valid for the interaction 
of particles in different tubes which is accurate when the transverse width, I, is much 
smaller than the intertube distance, d. Corrections to this picture have been discussed 
in references [21] and [25]. Below we will return to the question of finite transverse 
width when we treat bound states with more than one particle per tube. From here on 
we will adopt units and d for all energies and lengths. 

The potential above has the interesting property that for the case of n > 1 



and as a function of the angles we thus see that we can obtain both positive and negative 
values of the integrated interaction. The criterion for the existence of a two-body bound 
dimer in the limit of small U is a negative integral [271 EE] • We thus see that one can 
control the presence of a weakly-bound dimer by changing and 9. Below we will 
mostly address the situation where <p = tt/2 and 8 = 0, i.e. dipoles that are oriented 
perpendicular to the tubes. In this case, the integral in (j2]) is — ^ and this means 
that the system can produce bound dimers between any two particles that are located 
in different tubes for any n > 0. However, if a given few-body state has more than 
one particle in a single tube, then the interaction will be either attractive or repulsive 
depending on the sign of 1 — 3cos 2 (6 l ) cos 2 (0) (see (II]) with n = 0). Here we will only 
consider the regime where this quantity is positive, i.e. the case for which two dipoles in 
a single tube repel each other in order to avoid any collapsing states within the tubes. 

3. Few-body bound states 

We now present our results for up to five dipolar bosonic particles in multi-tube 
geometries. The discussion will involve analytical tools for addressing the limits of 
weak and strong dipolar interactions which will be compared to numerics using the 
exact stochastic variational approach [291 ED] • This allows us to determine the range of 
validity of the analytical methods that we employ. Along the way we will make a detailed 
comparison between the one- dimensional case with tubes and the two-dimensional case 
of multiple layers [7J [311 [321 [331 EH [351 ESI E7J. We will also compare to the results 
of McGuire [38] for the exact ground state energy of an iV-boson system with pairwise 
zero-range interactions in one dimension. Finally, we will investigate few-body states 
with more than one particle in a single tube for both the case of perpendicular dipoles 
(<f) = tt/2 and 8 = 0) and with tilted angles (8 = and <fi < tt/2). 
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3.1. Two-body states of two dipoles in two tubes 

The first case is one dipolar particle in each of two adjacent tubes with dipole moments 
oriented perpendicular to the tubes (0 = ir/2 and 6 = 0). This configuration has the 
Schrodinger equation 

, d 2 x 2 - 2 s 

U , 2 , 1W2 )$ = eS, (3) 



V dx 2 (X 2 + l) 5 /2 

where x is the relative distance of the two dipoles along the tube axis, e is the eigenenergy 
in units of and <3> is the wave function. We start with the analytically accesible limits 
of strong e/U ~ 1 and weak e/U «C 1 binding. 

VFeaA; binding. For a bound state we write the solution of ([3]) in the form 

*(«, i) = e— - r dy Sinh(K[g " y] V (y)$(K, y), (4) 

where e = — n 2 , k > [39]. Notice that we are 'normalizing' the wave function at x = oo 
and the integrating our way to any other value of x. For x < the term e~ KX appears 
to diverge, but the second term will also contribute and the combination yields the 
correct answer. For a potential that fall off faster than x~ l at infinity, we now obtain 
the equation for the binding energy 
i r°° 

dye«W(y)$( K ,y). (5) 



2 

For weak binding this can be solved iteratively for small k, i.e. 

j poo j poo poo 

« = - o / ^(j/)dj/ + - / dyy(y) / dy'(y - y')V{y') + ... (6) 

Here we have used that $(ac, y) ~ 1 + o(k 2 ). Upon inserting the potential in (JTJ with 
n = 1 and carrying out the integrals we find 

K = f/ + £/ 2 ^ + o(f/ 2 ), (7) 
lb 

which implies that to leading order in U, e = —k 2 oc — U 2 . These results were 
initially discussed by Landau [27] and Simon [28] but the procedure we present here 
allows simpler access to the higher-order terms. We note that the first term in (JH]) 
corresponds to a delta-function potential ~9^p^(^) with g = 2U (where we have 
reintroduced explicit units for clarity). Numerically we find that this potential gives 
good approximation for the real potential with U < 0.2. This is a valuable conclusion 
for studies that depend on a delta-function representation of the dipolar interaction in 
one dimension in order to describe many-body physics |21j . 

A clear difference to the 2D case is seen at this point. The weak-binding result in 
ID is e oc -U 2 , whereas in 2D it is e oc - exp(-8/£/ 2 ) [2HJ E33 El EE] . This is consistent 
with our expectation that weakly-bound ID dipolar few-body bound states are generally 
more stable than similar systems in 2D, and any sort of external perturbations would 
have a less severe effect in the ID setup. 

We can also consider the case where U < corresponding to one of the dipoles in 
the two tubes pointing in opposite directions which should be possible to achieve using 
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Figure 2. Energy, E, in units of as a function of dimcnsionless dipolar strength, 
U. Comparison between exact numerical solution of the Schrodingcr equation (j3|) and 
analytical predictions from table [1] for ID and 2D. The inset is a magnification of the 
small U regime. Note that the huge difference at small U between ID and 2D due to 
the exponential decrease of the binding energy at small coupling strength in 2D. 



AC fields [ID]. In this case the 2D setup with two adjacent planes will always have 
a bound dimer [30], whereas in ID we find numerically that it takes a finite strength 
U < —4.98 to bind the two-body system. We can understand the latter if we consider 
k = in (JSJ). This transforms the condition for the appearance of the bound state. For 
small U we have U + U 2 ^ = 0, which has two solutions; one corresponds to bound 
states with U > and another for bound states with U < — — ~ —5.09. The latter is 

7T 

within 2 percent of the numerical value. However, notice that the U < case opens up 
the possibility of binding and unbinding the dimer by tuning U as originally envisioned 
in reference |31j . 

Strong binding ID For this case we may assume that our wave function is strongly 
localised near the origin, i.e. (x 2 ) /e — > 0, such that the probability to be in the classically 
forbidden region is small. Under these conditions we can use standard techniques to 
obtain the binding energy and the wave function. First we decompose the potential 
near the origin 

2 o 00 

U . % ~^, =-2U + 6Ux 2 -45Ux 4 /4 + uYa n x 2n . (8) 

Using the assumption of strong localization we see that up to terms 0(-j=) we can solve 
the harmonic oscillator problem for the first two terms and include the third term via 
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perturbation theory. This yields 

e = -2U + V6U- 45/32 + 0(-^=). (9) 

A comparison of the numerical results for all values of U to the weak and strong 
binding expansions is shown in figure [2j Both the limit of weak and strong binding are 
well reproduced by the different approximation schemes. We see differences between 
numerical and analytical results only around U = 1. For completeness we also 
summarize all the analytical results and their regimes of validity in table [1] 





ID 


^JVdr 


-2U 


small U > 


e = -U 2 - f f/ 3 + o(f/ 3 ), ([/ < 0.2) 


large U > 


e = -2U + v/6Z7 - 45/32 + O(^), (£/ > 4) 




2D 


^JVdr 





small \U\ > 


e = -4exp(-2 7 - ^ + # - If + o(f/)), (17 < 0.8) 


large U > 


e = -2C7 + V2AU - 15/4 + O(^), (U > 7) 



Table 1. Comparison between ID and 2D binding energies of two dipolcs in two 
different layers or tubes with perpendicularly oriented dipolc moments. The regimes 
of validity of the different analytical formulae are indicated in the parenthesis. 



3.2. N-body Chains 

We now proceed to discuss the case where we have iV-body chains that consist of one 
dipolar particle in each of iV adjacent layers or tubes. In the 2D case, these sorts of 
structures have generated a lot of recent interest since one expects highly non-trivial 
few- and many-body dynamics [3H ESJ HU [42], H3j HU HH 06]. There is a similar interest 
in the one-dimensional multi-tube configurations [To*l [T9~| |2~T | |2"2"]. The Hamiltonian for 
the iV-body ID system is 

N 2 

# = Ef^+E y ( a ^)> ( 10 ) 

i=l i>j 

where Xij is the relative distance of the ith and j'th particles along the tube direction 
and Uij is the integer that signifies how far apart the tubes are that hold the ith and 
jth particles (n = 1 for adjacent tubes and so on). For the case where the pair-wise 
potentials of equal strength and zero-range, this Hamiltonian was studied in the classical 
papers of Lieb and Liniger [471 HE]- McGuire has shown that for this delta- function case, 
the bound state energy, E^, of the iV-body problem can be written very elegantly as 
En/E 2 = N(N 2 — l)/6, where E 2 is the two-body bound state energy of the same 
delta-function potential [38] . 

Here we consider a system of iV dipoles placed in iV different tubes with the 
potential in (fTOl) given by (CQ). This is intrinsically a case with long-range interactions 
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which has also generated some classical studies. In the particular case where the 
interaction behaves as r -2 this is the Sutherland-Calogero model (IHl [50] . In the 
present case with dipoles we have a long-range tail that behaves as r -3 . We would 
like to build a connection to the zero-range studies for the dipolar systems in ID. For 
general polarization angles, however, the interactions are anisotropic and more difficult 
to compare to the isotropic zero-range models. In this section we therefore study the 
case of perpendicular polarization (0 = 7r/2 and = 0) where the anisotropy is absent. 



20 l. 3 particles (toy model) 

4 particles (toy model) 

5 particles (toy model) 

3 particles (real chain) 

4 particles (real chain) 

5 particles (real chain) 

15 



LU 



5 - 



I 1 1 1 1 1 

2 4 6 8 10 

u 

Figure 3. Numerically determined -/V-body binding energies relative to the two-body 
energy, E 2 , as function of coupling strength, U, for iV-body chains using both a toy 
model where all interactions have n = 1 in (|Tj) and the realistic case where n is different 
for different pairs of particles in the chain. 

As we have demonstrated above, for U < 0.2 we get universal behaviour in a sense 
that the potential JJ ^-V = — 2C/5(|) accurately describes our system. Now we want 
to establish how well this zero-range approximation works on the few-body level with 
N > 2. To do so we first consider a toy model where all N particles interact with the 
same interaction corresponding to n — 1 in (JTJ. This is precisely the case where it makes 
sense to compare to the analytical result E^jE-i = N(N 2 — l)/6. In this way we can 
get a feeling for when the zero-range approximation works. For realistic chains n is not 
always equal to one. This means that these are easier to describe using a zero-range 
approximation since particle pairs that are located several tubes apart in a chain have a 
potential that effectively correspond to a smaller value of U (as compared to the n = 1 
case). 

The numerically calculated chain binding energies relative to the two-body energy, 
En/E 2 , for both the toy model and the real dipolar chains in ID are presented in 
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figure [31 We observe the delta-function behaviour in all cases for U < 0.15, and the 
strong coupling behavior is slowly approached with increasing interaction strength. In 
the strong coupling limit, the leading term for the toy model is En / 'E% = jV( ^~~ 1 ) which 
is simply the number of pairs that all contribute an energy E 2 . For the realistic chain, 
the expression is instead 

N N~k 

E N /E 2 = J2J2p= NH N-1 - H N-n (11) 
fc=l i=l 

where = YllLi ^ s the harmonic number of order k. We will elaborate more on 
the strong-coupling limit below. 

We see a large difference between the toy model and the real chain in figure [3] in 
the weak-coupling limit U <C 1. This is an important conclusion of our work that has 
implications for many-body studies on dipolar systems that work with delta-function 
approximations for the dipolar interactions. The toy model is relevant for an equidistant 
triangular tube configuration as studied in reference [2T]. Our results for this N = 3 
case show a very flat profile for all values of U and from the point of view of energetics it 
is thus a reasonable approximation to use the delta-function Zj£-V(x) = —2U5(^). The 
N = 4 and 5 would be applicable to a tube configuration with four and five nearest- 
neighbours respectively (the latter would of course not be possible with standard crystal 
lattices). However, in these cases the variation of the energies with U is much more 
drastic, indicating that a zero-range description is only good for very small U. For the 
real chain systems, the curves are rather flat for all N. This is connected to the fact 
that the terms with n > 1 in (11 01) will give a contribution that is suppressed. Due to 
the flatness of the curves in figure El we can get a good approximation for the iV-body 
binding energy of the realistic chains by using an effective model for iV-body systems 
with *j£-V(x) = -2U ef f(N)8(z) where 

Ueff(N) = _ 1} [NHSL, - Hgl,] U. (12) 

The blessing of working with a zero-range interaction comes at the price of having an 
iV-dependent coupling. This ensures that the iV-body bound states present in this 
geometry are properly described. We will use this effective dipolar strength to study 
the corresponding many-body system below. 

Chains in the oscillator approximation As we have seen for small U, the zero-range 
approximation gives reliable results. However, we would also like to discuss the strong 
binding limit where U 3> 1 and provide a simple analytical procedure for obtaining the 
energies and wave functions. This can be done using exact harmonic models [5TJ [52] 
where the full two-body interaction (the dipolar interaction in this case) is represented by 
a harmonic oscillator with parameters that are carefully chosen to reproduce energetics 
and structure at the two-body level [15]. We expect this to be a good approximation 
for U ^> 1 since the potential in ([T]) has a deep pocket in this case. However, it turns 
out that this approximation can be very accurate even for moderate U in the case of 
2D multi-layered systems with chains [53] . 
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Figure 4. Energy, E, in units of as a function of dimcnsionless dipolar strength, 
U. Comparison between exact numerical calculations obtained from the stochastic 
variational approach and the oscillator approximation described in the text for the 
case of perpendicular dipoles, i.e. <j> = it/2 and 6 = for the case of N = 3, 4 and 5. 



In figure H] we show a comparison of the oscillator approximation to exact numerical 
calculations for the case of perpendicularly oriented dipoles (<f) = tt/2 and 6 = 0). Note 
that the horizontal axis begins at U = 2 (for U < 2 the oscillator does not reproduce the 
exact results well). We get a very good agreement for iV = 3, while for higher particle 
numbers there is a clearly visible deviation. This stems from issues with the two-body 
interactions coming from the outermost tubes where the oscillator approximation starts 
to fail since the contributions are not well- described by an oscillator (the effective U for 
these terms is small for any U value considered here). However, the overall agreement 
is still within a few percent. One way to improve the agreement for higher particle 
numbers could be to introduce an effective three-body force (more generally an iV-body 
force) in addition to the two-body potential terms. Such an approach is very common 
in nuclear physics and in effective field theory studies of universal three-body physics 
[52(55]. 

We now change the angle 4> to investigate the effects of anisotropy. Figure [5] shows 
the case with = m and 9 = 0, where m satisfies cos 2 m = |. Recalling the discussion 
after (JTJ) this implies that the potential between two dipoles in the same tube vanishes. 
First we observe that the overall energy scale becomes smaller, because the depth of 
the potential and net volume of the potential is smaller than for the perpendicular 
polarization case in figure SJ For N = 3 and N = 4 we again see very good agreement 
between the oscillator approximation and the numerics. The largest chain with N = 5 
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shows what appears to be an even better agreement than for the <fi = 7r/2 and = case. 
What is also striking is that the oscillator results and the numerics cross around U ~ 5. 
This happens since our method ensures that the oscillator reproduces the slope of the 
numerical results in the limit U ^> 1 (as discussed for the 2D case in reference |53j). 
The optimization of the parameters for the two-body oscillator potential for 2 < U < 5, 
however, produces a slightly different slope and a crossing thus takes place. 
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Figure 5. Energy, E, in units of ' n jk- as a function of dimcnsionless dipolar strength, 
U. Same as in figure H] but for the case where (f> = <j) m and 8 = 0. This corresponds 
to vanishing intratube interaction, i.e. two dipolar particles in the same tube do not 
interact. 



3.3. Non-chain bound complexes 

We now consider the case where there are two particles in one tube. The simplest 
configuration of this sort is two tubes with three particles in total, i.e. one particle in 
one tube and two particles in the adjacent tube. At this point we need to take the 
width of the tube that we discussed in section [2] explicitly into account in order to work 
with an intratube interaction that best describes the experimental situation. This can 
be done by folding the potential with a gaussian wave packet of width / that represents 
the wave function of the system in the transverse direction (along the y-axis in figure [1]) 
This yields the intratube interaction 

(13) 

where 



/(*) 



: U rep X 3 (l — 3 cos 2 0cos 2 9)f(Xx), 
-2\z\ + v^tt(1 + z 2 ) exp(z 2 /2)erfc(|z|/v / 2) 



(14) 
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and A = d/l and U rep is the strength of the intratube interaction. In a system where 
the dipole moment is the same for all particles in all tubes we have U rep = U, where 
U is the strength of the intertube interaction. The dipole orientation and the induced 
dipole moment is typically controlled by applying an external electric or magnetic field 
that is constant across the whole system. However, if one applies also a field gradient 
then it should be possible to obtain a system where the induced dipoles moment varies 
from tube to tube (or layer to layer in the 2D multilayer case). In this case we can 
have a situation where U rep 7^ U and we thus vary both quantities independently in the 
current study. As discussed previously, we will consider only the repulsive case where 
1—3 cos 2 <p cos 2 9 > in order to avoid possible intratube collapse of the system. Unless 
explicitly states, we use A = 5 everywhere (see reference [25] for a discussion of the 
effects of changing A). 

We start by considering the case where all particles are perpendicularly polarized, 
i.e. = 7r/2 and 9 = 0. Observe that if U rep = we have a trimer energy that is below 
the dimer energy, and if U rep — > 00 we can not bind three particles. We now want to 
find the value at which the three-body system is bound for a given value of U, so 
we keep U fixed and vary U rep . In reference (36] it has been shown that U^ p /U is an 
increasing function of U by an argument that does not depend on the dimensionality 
and thus applies equally well to layers in 2D and tubes in ID. In turn we only need to 
consider U — > 00 to obtain an upper bound for U^ p /U. This upper bound is a decreasing 
function of A and for A = 5 we get U^ r ep /U < 0.41. In order to find the lower bound for 
Url p /U we will use the variational principle with a normalised three-body wave function 
of the form ty(x,y) = <f>(x)(p(y), where is the normalised dimer wave function and 
we have used Jacobi coordinates \[2x = x\ — £2 and a/3/2?/ = x 3 — {x\ + x%)j1. From 
this we obtain the upper bound for the three-body energy 

El ar = (q\H 3 \^)=E 2 + 



^|(r w + J ^(^[^(l^-x+^yD + niy-x-^/^Djda;)^)^^ 

where the kinetic term for the y coordinate is T y = —^-^p- Here V rep is given in 
(I13p and V is given by ([I]) with n = 1. Equation (jSJ) in principle allows us to estimate 
where E^ ar < E2 holds, which gives us the lower bound for U^ r ep . However, in general this 
equation is not easy to analyse so we obtain the lower bound using the simpler condition 
for the effective potential acting on ip in (I15p which is that it must have negative net 
volume. The limit of zero net volume requires 

$ 2 (*) \Vre P {\\J\x + + V{\^-x - ^y\)]dxdy = 0, (16) 

which can be simplified to 

POO 

V rep (y)+V(y) dy = -2U+U rep X 2 / f(y)dy = -2U+X 2 U rep = 0.(17) 




This lower bound is a decreasing function of A, and for A = 5 we get U^ p /U > 2/25. 
Collecting our results, we have for A = 5 that 0.41 > U° l p /U > 2/25. This means 
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that if we could manipulate the values of the intratube and intertube interactions 
independently, we could create or destroy the three-body bound state in two tubes 
with perpendicular polarization. This effect is completely due to the finite range of 
the potential, because for zero-range potentials any value U rep /U will produce a bound 
state (the system becomes equivalent to the fermionic states studied by McQuire in 
reference [56] for which a bound state always exists). 

0.34 | 1 1 1 1 1 




0.18 1 1 1 1 1 1 

10 20 30 40 50 

u 



Figure 6. Ratio U^ p /U as a function of U at which a three-body state with two 
dipolcs in one tube (plane) and one in the adjacent tube (plane) becomes possible 
in ID (2D). The three-body state exists below the lines in the plot. The lines drops 
sharply to zero for small U. 

In figure El we show the critical value, U^ p /U, as a function of U for both the case of 
ID tubes and for comparison we also plot the results of the same study for a bilayer in 
2D [36]. Both of these calculations have been done with A = 5. We see a clear difference 
between the critical values for three-body bound state formation in the two cases with 
the ID situation allowing for larger repulsive interactions. This is connected to the price 
one pays for localizing particles in a bound states complex. In 2D we need to localize 
in two independent directions which is more expensive than ID. As expected, ID tubes 
are more favourable for formation of these larger complexes. 

Having discussed how a tuning of the intratube repulsion can bind a three-body 
state, we now consider an alternative route to bind this state which is by tilting the 
angle, 0. We assume that U rep = U, and are looking for values of <fi and U for which a 
bound state exists. For a given angle <fi we introduce the notation U cr (<p) which is the 
value at which the two-body dimer (one particle in each of two adjacent tubes) energy 
is equal to the three-body energy which defines the threshold. We start with three 
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particles as before but tilted by angles 9 = 0, G [0 m , 7r/2]. We already proved that for 
</> = 7r/2 we do not have a bound state if U rep = U . However, we know that for = m 
we always have a three-body bound state, because the intralayer interaction vanishes. 
This implies that there is a range of angles where the three-body state exists. 

To approach this problem analytically, we can calculate the net volume as above. 
This yields 

dy = 2[/(cos 2 - 1) + U(l - 3 cos 2 0)A 2 . (18) 

Putting it to zero we obtain the smallest angle below which the three-body bound 
state exist for any value of U . The condition for a three-body bound state is 
< cos-^W^^) which in the case of A = 5 yields < 0.975 or 55.9°. In the strict ID 

limit where A — > oo, we find the result < m . The function g^" 2 is decreasing with A 
which is a result of the behaviour of the function in ( fl4l) which has a finite value for zero 
argument for any finite A. This means that the intratube potential has been regularized 
by the finite width of the transverse confinement and does not have the strict ID form 
x~ 3 expect in the limit of A — > oo. In strict ID we thus expect that the repulsion will not 
allow this three-body bound state due to the cubic singularity. For the more realistic 
case with a finite A we can have bound states for angle in the regime > m as well. 
An additional bound can be found by using classical arguments, i.e. neglecting kinetic 
energy in the Schrodinger equation and studying just the potential landscape [36]. This 
yields a bound of > 1.15 or 66° so that for any U we do not have three-body bound 
states in the case where U rep = U. In summary, we have found that for m < < 1.15 
there exist a value of U cr (<j)) G [0, oo), such that for U > U cr (4>) the bi-tube trimer is 
bound. 

In figure [7] we plot U cr as a function of for three particles in two tubes in the case 
where U rep = U, i.e. the case of current experimental interest where the dipole strength 
is controlled by one external field that is constant over the entire samples. The points 
indicate the values at which the dissociation points have been obtained numerically 
(connecting lines are a guide to the eye). The three-body bound states are present on 
the left-hand side in the diagram since here the intratube repulsion is not strong enough 
to break the bound state apart. The results obtained here are consistent with results 
reference [25] for U < 15 where a comparison can be made. For U > 15 we find a very 
sharp increase in U cr (<f)) that takes place at ~ 1.1, a few degrees before a classical 
calculation renders the system unbound at = 1.15. 

In order to address the question of additional tubes and their effect on bound states 
with more than one particle per tube, we have also consider a four-body state in three 
tubes (as in figured]) with two dipoles in the middle tube and one in each of the two outer 
tubes. The results are also presented in figure UJ We again see a sharp increase for large 
before we get to the classically unbound value of = 1.34 or 76.5°. However, there is 
a significant increase of region in which bound state with two particles in a single tube 
are present due to the extra attraction from adding a particle in an adjacent tube on the 
other side of the complex as compared to the three-body system with two and one in two 
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Figure 7. Critical value, U cr (4>), where the trimer in two tubes (one and two particles) 
and tetramer in three tubes (one, two and one particles) dissociates as function of the 
tilting angle 4> for 9 = 0. The dissociation thresholds where calculate at the dots 
and squares with connecting lines to guide the eye. The vertical lines mark the angles 
beyond which a classical calculation (neglecting kinetic energy for all particles) predicts 
no bound states can occur because the system is overall repulsive. 

tubes. Other four-body states can occur in two tubes as discussed in reference [25], but 
the parameter regime is severely suppressed due to the extra repulsion in the three-plus- 
one and two-plus-two four-body states in two tubes. We therefore expect that the case 
of three tubes could have interesting 'quasi'- Wigner crystal states for large U beyond 
the clustered states studied in reference [18] where the clusters are not of the same size 
in each tube (here with cluster size one in the outer tubes and two in the inner tube). 
We can continue this kind of game for five tubes in a one-two-one-two-one configuration 
and so on. 

4. Scattering 

While we are mainly concerned with negative energy bound state structures in this 
paper, we now present a short overview of some scattering results that can be deduced 
directly based on the formalism we introduced above. One basically only needs to make 
an analytical continuation of the bound state discussions above in order to arrive at 
results for scattering of dipolar bosons in ID tubes. Since current experiments with 
cold molecules aim for the low temperature regime where quantum degeneracy can be 
reached, we will focus exclusively on the low-energy scattering regime here. 
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The formal solution of the scattering problem, satisfies the equation 

= e ikx ~ fdy S]n{k[X k ~ v]) V(yMk : y) : (19) 

where the boundary condition is an out-going flux to the right with coefficient one of the 
form e tkx [39]. Here we define the reflection, R(k), and transmission, T(k), coefficients 
through ty(k,x — > — oo) = T(k)~ 1 e tkx + ^|^e _ifcx to be consistent with standard ID 
potential scattering. By taking the limit x — >■ — oo we obtain 

2ik 

T ^ ~ 2%k - dxV(x)e- ikx ^(k, x) ' ^ 

For k = in, k > 0, poles of T{k) define bound states with the binding energy defined 
through the equation 

^e) = n = --J dye^V^^y), (21) 

which is consistent with ([5]) and allows us to use the same techniques as before to obtain 
iterative improvements of the solution. 

Consider now the case of dipoles oriented perpendicular to the layers, i.e. (j) = n/2 
and 6 = 0. In the Limit of weak binding and low-energy scattering, [/< 1 and kd <C U, 
we find for the transmission coefficient of two particles in the two different tubes 

T = i-, (22) 

and for the case of two particles in the same tube we find 

kd 

These results are deceptively simple, yet they contain all the information about the 
dipolar physics in one- dimensional tubes at low energy in the small U limit. As expected 
they vanish with k as k — > which is a manifestation of the Wigner threshold law. Also, 
as Kb oc U in the weak-coupling limit (see table [1]), we see that both expressions become 
large at very small U. This makes sense since the potential becomes very weak (and 
likewise the derivative of the potential). We also see that in the strict ID limit, A — > oo, 
the transmission is completely suppressed and total reflect is expected. This is the 
emergence of the impenetrable boson regime called the Tonks-Girardeau limit j57J EE] • 
This limit can be studied with dipolar bosons in a ID setup with an (in-tube) harmonic 
trap as discussed in reference [59] . 

Since we would like to obtain an effective model for dipolar systems of iV-chains 
using zero-range interactions, we need to relate the transmissison coefficients above to 
its zero-range equivalent. If we write the potential as V(x) = gS(x), we find 

2ik 2ikh 2 , , 

r m = ^-^ - -— , (24) 

for low-energy scattering. We see that we recover the g = —2U result for the case of 
two dipoles in different tubes (using our units of energy and length given in section [2]) 
and we find g = UX 2 for dipoles in the same tube. 



T = -ijTVf (23) 
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5. Many-Body Physics 

We now want to address the case of a many-body system of N adjacent tubes with M 
particles in each tube. We will focus on the case of perpendicular orientation of the 
dipoles with respect to the tubes (0 = n/2 and 6 = 0) and comment on the case of 
general angles in the end. In the perpendicular case, we expect no bound complexes 
beyond the chain structures that are bound states of N bosons with one in each of N 
adjacent tubes. Recent Monte Carlo studies of the many-body system indicate that the 
longest chains possible, the iV-chain for N tubes, are the relevant degrees of freedom in 
these systems [221 IH] and we therefore focus our attention on these chains. 

The results that we have presented thus far show that one can make an effective 
model of chains using zero-range interactions with suitably adjusted coupling strength 
as in (112j) . This is extremely useful for many-body physics since a number of ID N- 
body models with zero-range interactions are exactly solvable as discussed above. For 
the case of bosons with attractive zero-range interactions one in fact finds a gas-like 
state similar to the Tonks- Girardeau gas of repulsively interacting bosons, denoted 
the super- Tonks- Girardeau gas [60l I6T1 [62] . For both attractive and repulsive zero- 
range interactions these systems are also possible to describe using bosonic Luttinger 
liquids (see reference [63] for a recent review of ID bosonic systems). These studies 
have been boosted by recent experimental success in realizing these ID quantum gas 
[M| [65] . For dipolar bosons this Luttinger liquid behavior was studied by several authors 

Here we want to address the many-body behaviour by using the few-body physics 
information that we have derived above. For N adjacent tubes with M particles in each, 
we can build a system of M chains, and then we may study the effective inter-chain 
interactions. Using a zero-range interaction model allows us to map the problem into 
the realm of exactly solvable ID systems and/or Luttinger liquid studies and thus infer 
the many-body dynamics. We will stick to the low-energy phases here in order to use 
the low-energy scattering and weak-coupling results we have derived above. 

The effective zero-range coupling in (JT2~]) has the property that it (approximately) 
reproduces the iV-body binding energy of a single chain. However, if one considers 
two chains then one realizes that this effective interactions also describes the attractive 
force between two different chains. This follows readily from the fact that the intertube 
interactions that are responsible for the chain-chain interactions are the same as those 
that provide the intertube attraction in a single chain. If we ignore for a moment any 
intratube repulsion, we can estimate the ground state energy of a system with M chains 
of length N to be simply 

E M = M{M2 ~ 1) E CC (N) + ME CC (N), (25) 
6 

where E CC (N) is the two-body energy between two different iV-chains (consider as 
effective degrees of freedom in the system) calculated using the potential ^V a (x) = 
—2U e ff5{^). The first term is obtained by mapping to the result obtained by McGuire 



Bound states of Dipolar Bosons in One- dimensional Systems 



18 



for M bosonic chains and the second term comes from the internal energy of each 
of the M chains of length N. The iV-dependence is given implicitly through U e ff 

The repulsive intratube interaction has thus far been ignored. However, it is 
not difficult to take into account through a zero-range interaction term of the form 
1! ^-V r {x) = iV£/A 2 5(|). The factor of N is obtained by noting that there will be one 
repulsive contribution from each tube seperately. In our effective model of chains, we 
thus have a competition between two interactions of opposite sign, V a (x) and V r (x). We 
can use this to esimate when the system of chains constitute an attractive or a repulsive 
Bose gas by determining when V a (x) = V r (x). Since both are proportional to U, this 
will give a purely geometric condition. If we define 

HN, A) = N{ ™_ 1} [NHP_ 1 - ff« J - NX 2 , (26) 

then V a (x) = V r (x) corresponds to F — 0, while F > gives an attractive and F < 
a repulsive chain system. We thus see immediately that except for very small A, the 
system will be in the repulsive regime for any N > 2. Small A corresponds to a soft 
transverse confinement that requires careful treatment beyond our quasi-lD formalism. 
We therefore conclude that the system of dipolar chains in the weak-coupling limit 
with dipoles oriented perpendicular to the tubes should behave as a repulsive Tonks- 
Girardeau gas. The effective parameters introduced here can be used to estimate ground 
state properties using the Lieb-Liniger model [4"Tll4"8"] and/or Luttinger liquid formalism 
can be used to obtain correlation functions etc. 

If we tilt the dipolar orientation away from perpendicular, we expect a smaller 
effective coupling due to the angular factor in (JTJ). The repulsive part will, however, 
remain dominant and we expect the same conclusion as above to hold for small tilting. 
If we go to the value <fi m and 9 = where the intratube repulsion is zero, we naturally 
get an attractively interacting system. However, here we need to worry about more 
complicated complexes than bound states of chains as we have discussed above (this 
was first pointed out in reference [24]). One could imagine that if we take a long 
chain (large N), then it may be possible to attract one more particle to bind to this 
chain around the middle tube. This may open the possibility for scattering events 
between chains where this extra particle is exchanged, akin to transfer reaction known 
for instance from nuclear physics. The nature of potentially more complicated bound 
structures constitute an interesting direction for future work. 



6. Conclusions and outlook 



We have studied dipolar bosonic particles confined to a setup consisting of a number of 
one-dimensional tubes. For the case of equidistant adjacent tubes, we calculated the few- 
body bound state structures for up to five particles using both analytical and numerical 
tools in the weak and strong binding limits. For the case of perpendicularly oriented 
dipoles, we find that chain with one particle in each of a number of adjacent tubes are 
the most stable structures, and that more complicated bound structures with multiple 
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particles in a single tube are unlikely unless one tilts the dipole orientation away from 
the perpendicular direction. The bound states in our study should be directly observable 
using experimental techniques such as lattice shaking [66], RF spectroscopy [671 EE] or 
in-situ optical detection [2^ [25]. 

Our results for the weak-coupling limit can be used to calculate the low-energy 
scattering properties of the system and make an effective zero-range model for bound 
chains of dipolar bosons that takes both attractive and repulsive terms into account. 
This allows us to use our knowledge of the few-body physics to address the many-body 
physics of these systems which can be mapped onto exactly solvable model in ID and 
Luttinger liquid dynamics. Estimation of the magnitudes of the competing attractive 
and repulsive terms in a system of chains indicate that for perpendicular orientation 
the repulsion is always dominant and a Tonk- Girardeau type of system should emerge. 
This could be changed by tilting the angle of the dipoles with respect to the tubes 
which changes the balance between repulsive intratube and attractive intertube terms. 
However, the potential presence of bound state complexes with two or more particles 
per tube need to be considered in this case. The latter will be the focus of future 
investigations. 

In the limit of strong binding we demonstrated that models based on replacing the 
true dipolar potential by a harmonic interaction are quite accurate for describing the 
energetics. This yields exactly solvable iV-body models that could for instance be used 
to study thermodynamic properties in ID dipolar systems [691 ED] • Another intriguing 
possiblity is to study system where one or more of the tubes have perpendicular dipoles 
that point in the opposite direction. These flips can be obtained by using AC fields [71] . 
For two dipoles in two adjacent tubes with opposite orientation of the dipoles we proved 
here that in the limit of weak dipole moment, there will be no bound two-body state. 
However, this does no in general rule out the possibility of a bound three-body state, 
a so-called Borromean state [72]. We could thus imagine a regime where three-body 
states are the lowest non-trivial bound states and thus give rise to an effective gas of 
trimers. While all of the questions addressed in this paper have been concerned with 
bosonic dipolar particles, we expect that similar bound state structures should arise 
when using fermionic dipoles. The fermionic case should correspond closely to studies 
on one-dimensional multi-component systems [731 E] if one maps the component index 
onto tube index. 
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